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ABSTRACT 

Aims. We propose an explanation for the considerable scatter of the abundances of neutron capture elements observed in low-metallicity stars 
in the solar vicinity, compared to the small star- to- star scatter observed for the Qf-elements. 

Methods. We have developed a stochastic chemical evolution model in which the main assumption is a random formation of new stars subject 
to the condition that the cumulative mass distribution follows a given initial mass function. 

Results. With our model, we are able to reproduce the different spreads of neutron capture elements and Qf-elements in low-metallicity stars. 
Conclusions. The reason for different observed spread in neutron capture elements and Qf-elements resides in the random birth of stars, coupled 
with different stellar mass ranges, from which Qf-elements and neutron capture elements originate. In particular, the site of production of a- 
elements is the whole range of massive stars, from 10 to 80M o whereas the mass range of production for neutron capture elements lies between 
12 and 30M o . 
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1. Introduction 

Early work by Gilroy et al. (1988) first proposed that the stel- 
lar abundances of very heavy elements with respect to iron, 
particularly [Eu/Fe], show a large scatter at low metallicities. 
Their work suggested that this scatter appears to diminish with 
increasing metallicity. This was confirmed by the large spread 
observed later on in the [Ba/Fe] and [Eu/Fe] ratios in halo stars 
(e.g., McWilliam et al. 1995; Ryan et al. 1996). A more ex- 
tensive study by Burris et al. (2000) confirmed the very large 
star- to- star scatter in the early Galaxy, while studies of stars at 
higher metallicities, involving mostly disk stars (Edvardsson et 
al. 1993; Woolf, Tomkin, & Lambert 1995), show little scat- 
ter. In the last few years, a great deal of observational work 
on galactic stars appeared: Fulbright (2000); Mashonkina & 
Gehren (2000, 2001); Koch & Edvardsson (2002); Honda et al. 
(2004); and Ishimaru et al. (2004). All these works confirmed 
the presence of the spread for these elements and this spread 
can reach 2 dex at [Fe/H^l — 3. It is worth noting that such a 
high spread is not found for the [a/Fe] ratios in very metal-poor 
stars (down to [Fe/H] = -4.0, Cayrel et al. 2004); in fact, the 
spread for these elements is around 0.5 dex at low metallicities. 
This fact suggests that the spread, if real, is a characteristic of 
these heavy elements and not only due to an inhomogeneous 
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1 We adopt the usual spectroscopic notations that [A/B] = 
log m (NA/N B )* - logi (N A /N B ) Q and that loge(A) = log m (N A /N H ) + 
12.0, for elements A and B 



mixing in the early halo phases, as suggested by several au- 
thors (Tsujimoto et al. 1999; Ishimaru & Wanajo 1999). 

Recently, several studies have attempted to follow the en- 
richment history of the Galactic halo with special emphasis 
given to the gas dynamical processes occurring in the early 
Galaxy: Tsujimoto, Shigeyama, & Yoshii (1999) provided an 
explanation for the spread of Eu observed in the oldest halo 
stars in the context of a model of supernova-induced star forma- 
tion; Ikuta & Arimoto (1999) and McWilliam & Searle (1999) 
studied the metal enrichment of the Galactic halo with the help 
of a stochastic model aimed at reproducing the observed Sr 
abundances; Raiteri et al. (1999) followed the Galactic evolu- 
tion of B a by means of a hydrodynamical N-body/smoothed 
particle hydrodynamics code; and Argast et al. (2000) concen- 
trated on the effects of local inhomogeneities in the Galactic 
halo produced by individual supernova (SN) events, account- 
ing in this way for the observed scatter of some (but not all) el- 
ements typically produced by type II SNe. Besides a spread for 
r-process elements, they also predicted a spread of more than 2 
dex for Mg and O at [Fe/H] =-3, which is too large compared 
to the observational data. Finally, Travaglio et al. (2004) inves- 
tigated whether incomplete mixing of gas in the Galactic halo 
can lead to local chemical inhomogeneities in the interstellar 
medium (ISM) involving heavy elements, in particular Eu, Ba, 
and Sr. They reproduced the spread for these elements, but left 
unexplained why the spread is present only for neutron capture 
elements whereas it is much smaller for the other elements (for 
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example ^-elements). Concerning the ^-elements, Argast et al. 
(2002) and Karlsson & Gustafsson (2005) attempted to explain 
the small dispersion in these elements from two different points 
of view: a different iron yield distribution as a function of pro- 
genitor mass for Argast et al. (2002), and cosmic selection ef- 
fects favoring contributions from supernovae in a certain mass 
range for Karlsson & Gustafsson (2005). 

The purpose of this paper is to explain the different spread 
observed in the halo stars for different elements. Cescutti et al. 
(2006) have suggested the best stellar yields to fit the trend of 
the mean abundances for Eu and Ba in the framework of a ho- 
mogeneous chemical evolution model (Chiappini et al. 1997, 
2001). In the present paper, using the same yields of Cescutti 
et al. (2006), we show the results of a stochastic chemical evo- 
lution model that we have developed to explain the different 
observed spreads for different elements in the halo stars. The 
paper is organized as follows: in Sect. 2 we present the obser- 
vational data; in Sect. 3 the inhomogeneous chemical evolution 
model is presented; and in Sect. 4 the adopted nucleosynthesis 
prescriptions are described. In Sect. 5 we present the results; 
and in Sect. 6 some conclusions are drawn. 

2. Observational data 

For the extremely metal-poor stars (-4<[Fe/H]< -3), we 
adopted the recent results from UVES Large Program "First 
Star" (Cayrel et al. 2004, Francois et al. 2007). For the abun- 
dances in the remaining range of [Fe/H], we took high quality 
data of various sources published in the literature. These are 
are sumarized in Table [T] We homogenized all of these data 
by normalizing them to the solar abundances by Asplund et al. 
(2005). 

3. Inhomogeneous chemical evolution model for 
the Milky Way halo 

For comparison with a homogeneous chemical evolution 
model, we choose to use the same parameters of the chemi- 
cal evolution (star formation rate, initial mass function, stellar 
lifetime, nucleosynthesis, gas density threshold) as those of the 
homogeneous model used in Cescutti et al. (2006). For this rea- 
son, we model the chemical evolution of the halo of the Milky 
Way for the duration of 1 Gyr. We consider that the halo con- 
sists of many independent cubic regions, each with the same 
typical volume, and each region does not interact with the oth- 
ers. We decided to have a typical volume of 2.8 • 10 6 pc 3 , in this 
way the surface of the volume, taken as the side of a cube, is 
2 • 10 4 pc 2 . The number of assumed volumes is 100. In order to 
ensure good statistical results, we test that a larger number of 
volumes and foud that the statistics do not significantly change, 
but the model is slowed down. The dimension of the volume is 
large enough to allow us to neglect the interaction among differ- 
ent volumes, at least as a first approximation. In fact, for typical 
ISM densities, a supernova remnant becomes indistinguishable 
from the ISM - i.e., merges with the ISM - before reaching 
~ 50pc (Thornton et al. 1998), less than the size of our cubic 
region. We do not use larger volumes because we would lose 
the stochasticity we are looking for; in fact, as we tested for 



increasing volumes, the model tends to be homogeneous. We 
use timesteps of lMyr, which is an interval of time shorter than 
any stellar lifetime considered here; in fact, the maximum stel- 
lar mass considered here is 80M© with a lifetime of ~3Myr. At 
the same time, this timestep is longer than the cooling time of 
the SN bubbles, which is normally ~ 0.1-0.2Myr and at maxi- 
mum 0.8 Myr. We describe later in this same Sect, how we cal- 
culate this cooling time. With this timestep we can say, in first 
approximation, that an homogenous mixing takes place and, at 
least for the first period, only a few stars pollute the ISM; this 
is again important from the point of view of the stochasticity of 
the model. 

In each region, we assume the same law for the infall of 
the gas with primordial composition; originally, in the homo- 
geneous model, we have this law for the halo phase: 



dt 



nfalr 



-t/T 



(1) 



where the parameter I^ all is equal to \6M Q Gyr~ l pc~ 2 , and this 
law is a function of surface gas density. However, we want to 
use this law for our box, and so considering the side of the cube 
(2 • 10 4 pc 2 ) as the typical surface we obtain in each volume, 
multiplying the previous equation by this typical surface: 



U1VI tOt W _ Vol -t/T 
dt ' n f aU 



(2) 



where ff aU is 320M o Myr- 1 ; r is in both cases lGyr. We con- 
sider a threshold in the gas surface density, which means that 
below this threshold the star formation is blocked, due to the 
insufficient gas density in the system. The value of this gas sur- 
face density threshold is cr thresho i d = 4M Q pc~ 2 , as in the homo- 
geneous model, so with the surface we consider, it becomes in 
each volume M t h re shoid = 80 • 10 3 M o . The star formation rate 
if/(t) is defined: 



m = y( 



M gas (t) 
^threshold 



1.5 



(3) 



where M gas (t) is the gas mass inside the considered box 
and the parameter v is 80M©Mvr _1 ; this parameter is set to 
match the star formation rate of the homogenous model during 
the halo phase. So, when the gas density threshold is reached 
for the first time, the mass that is transformed into stars in a 
timestep (hereafter, M n ™ rs ) is 80 M Q . In this way, we note that 
the code can form at any timestep stars of any mass up to the 
considered maximum stellar mass (i.e. 80 M©), and the code 
can form all the stellar masses up to 80 M©. So, with the con- 
sidered surface (with the considered gas density threshold), we 
are also overcoming the risk of having a bias toward low mass 
stars, when the code randomly chooses the masses of the form- 
ing stars. 

Knowing the M n ™ rs , we assign the mass to one star with a 
random function, weighted according to the initial mass func- 
tion (IMF) of Scalo (1986) in the range between 0. 1 and 80 M Q . 
Then we extract the mass of another star and we repeat this cy- 
cle until the total mass of newly-formed stars exceeds M n ™ rs - 
In this way, in each region, at each timestep, the M n ™ rs is the 
same but the total number and mass distribution of the stars 
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are different, and we know the mass of each star contained in 
each region, when it is born and when it will die, assuming 
the stellar lifetimes of Maeder & Meynet (1989). We compute 
the chemical evolution in the following way: at the end of its 
lifetime, each star enriches the interstellar medium with newly- 
produced elements (see the next Sect.) as a function of its mass 
and metallicity. The total mass of each element is determined 
at the end of the lifetime of each star, taking into account the 
enrichment due to that star and due to the already present mass 
of each element locked in that star when it is born. The model 
does not take into account the pollution produced by stars with 
mass < 3M© because their lifetimes exceed the duration of the 
simulation. The existence of SNe la is also taken into account, 
according to the prescriptions of Matteucci & Greggio (1986), 
in the single degenerate scenario. We consider that a fraction of 
the stellar masses in the range 3-16M© would be binary systems 
which can produce SN la. The correct number of observed SN 
la at the present time (see Cappellaro et al. 1999) is reproduced 
if this percentage is 5% and we evaluate, again using a random 
extraction, if the system is a progenitor of a SN la. We also 
determine the mass of the secondary through the distribution 
function: 

/0u) = 2 1+ ^(l+ r )^ (4) 

with y = 2, where [i is the fraction of the mass of the secondary 
star to the total mass of the system and in this way we know 
when the SN la will explode; in fact, its lifetime is that of the 
secondary star of the binary system. In Fig. [T] the SNe la rate 



is compared to the SNe II rate. Due to the threshold in the gas 
density that we impose, the star formation starts only after 250 
Myr and, at this stage, the first SNe II start to explode. With a 
time delay of about 30Myr, the first SNe la take place. In our 
model we consider the single degenerate scenario and 30Myr is 
the shortest timescale for a SN la to explode, being the lifetime 
of the most massive progenitor leading to the formation of a C- 
O white dwarf (i.e. a star of 8M Q ). We stop the star formation 
at 1 Gyr and the SNe II rate falls abruptly to zero, correctly, 
the largest lifetime of a SNII being 30Myr, whereas the SNe la 
continue to explode, the lifetime of the progenitors of a SN la 
being as long as 10 Gyr. 

We assume that the ejecta produced by the dying stars mix 
with the surrounding ISM only after these ejecta cool down 
significantly. The cooling timescale of the ejecta is calculated 
as follows: 

_ 3kT 

cooling ~ nA(T,z) ( } 

where n, T, and z are the average number density, the tem- 
perature, and the metallicity in the cubic region, respectively. 
A(r, z) is the cooling function, taken from Sutherland & Dopita 
(1993). We know the total gas mass and the thermal energy (as- 
suming 10 51 erg of energy after each SN explosion), therefore 
in each cubic region, at each timestep, we can calculate n and 
T. The resulting cooling timescales are normally ~ 0.1-0.2Myr 
and reach at maximum 0.8 Myr in a set of 100 volumes. Since 
the cooling timescales are always smaller than the timestep of 
our simulation (1 Myr), we decided to neglect this delay time, 
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time (Myr) 

Fig. 1. The SNIa rate (solid line) and SNII rate (dashed line) in 
the halo. Note that the rate of SNIa is multiplied by ten. 



so that the mixing in each simulated cubic volume is instan- 
taneous. The model follows the chemical evolution of 10 ele- 
ments (Si, Ca, Mg, Fe, Sr, Y, Zr, Ba, La and Eu) in each region. 
If the model is correct, our predictions will approximate the re- 
sults of the homogeneous model, having the same parameters, 
as the number of stars increases. On the other hand, our model 
shows the spread of chemical enrichment at low metallicity that 
can be produced by different stars of various masses, where the 
number of stars is low and the random effects in the birth of 
stellar masses are important. 

4. Nucleosynthesis Prescriptions 

For the nucleosynthesis prescriptions of the Fe, Si, Ca, and 
Mg, we adopted those suggested by Francois et al. (2004) both 
for single stars and SNela. They started with the Woosley & 
Weaver (1995) yields for massive stars at the solar composition, 
metallicity independent; they use the Iwamoto et al. (1999) 
yields for type la SNe. Then, they compared the model re- 
sults with observational data with the aim of imposing con- 
straints upon stellar yields. To best fit the data in the solar 
neighborhood, they modified the original yields by Woosley & 
Weaver (1995) and by Iwamoto et al. (1999). In particular, they 
changed Mg by increasing in stars with masses from 1 1 to 20 
M Q and decreasing in masses larger than 20 M . The Mg yield 
has also been increased in SNe la. The Si yields have increased 
slightly in stars above 40 M©. Ca and Fe (the solar abundance 
case) in massive stars from Woosley & Weaver (1995) are the 
best to fit the abundance patterns of these elements since they 
do not need any changes. We underline that the site of produc- 
tion of a-elements is the whole range of massive stars. For Fe, 
the main producers are SNela, but a fraction arises also from 



the whole range of massive stars. For the nucleosynthesis pre- 
scriptions of the r-process contribution, we used those of model 
1 for Ba and Eu (see Cescutti et al. 2006), and the results of 
Cescutti et al. (2007a) for La. For Sr, Y, and Zr we used the 
results shown in Cescutti (2007b). We have assumed that these 
neutron capture elements are produced by r-process in massive 
stars, but only in the mass range 12-30 M©. These empirical 
yields have been chosen to reproduce the surface abundances 
for all these neutron capture elements in low-metallicity stars 
as well as the Sr, Y, Zr, Ba, Eu, and La solar abundances, tak- 
ing in account the s-process contribution at high metallicities 
by means of the yields of Busso et al. (2001) for La and Ba and 
those of Travaglio et al. (2004) for Sr, Y, and Zr in the mass 
range 1.5-3M©. Eu is considered to be a purely r-process el- 
ement produced in the same range of masses. This choice is 
extensively discussed in Cescutti et al. (2006) and the nucle- 
osynthesis prescriptions of the r-process contribution are sum- 
marized in Table [2] In this model, we do not take into account 
the s-process contribution for these elements because its en- 
richment timescale exceeds the duration of the simulation, this 
process taking place in low mass stars (1.5-3M©). 

5. Results 

5.1. The ratios of a -elements and neutron capture 
elements to Fe 

We discuss here the results of our simulations compared to 
the observational data and to the prediction of the homoge- 
neous model of Cescutti et al. (2006). We show the [Eu/Fe], 
[Ba/Fe], [La/Fe], [Sr/Fe], [Y/Fe], and [Zr/Fe] ratios as func- 
tions of [Fe/H] in the Figs. [2 [SHOO and[7l respectively; for 
a-elements we show [Si/Fe], [Ca/Fe], and [Mg/Fe] as functions 
of [Fe/H] in Figs. [8j [9j and[T0l respectively. In these Figs., we 
show the simulated living stars at the present time (blue points) 
together with the observations (red triangles). 

It is worth noting that our model reproduces the large 
spread observed in the abundances of metal-poor stars for the 
neutron capture elements and, at the same time, the small 
spread for the a-elements. The chemical enrichment observed 
in very metal-poor stars is due only to nucleosynthesis in mas- 
sive stars. The site of production of the a-elements is differ- 
ent than the one of neutron capture elements (as described in 
the previous Sect.): the a-elements and Fe are produced in the 
whole range of massive stars (10-80M©); on the other hand, 
the neutron capture elements are assumed to be produced in 
the smaller range between 12 and 30 M©. Therefore, in re- 
gions biased randomly toward stars less massive than 30 M©, 
the ratio of neutron capture elements over Fe is high. The op- 
posite happens in regions where a large fraction of the stars are 
more massive than 30 M©. This fact produces, in our inhomo- 
geneous model, a large spread for the rates of neutron capture 
elements to Fe, but not for the ratio of the a-elements to Fe, 
since the a-elements and Fe are produced in the same range. 
The idea of a small range of the r-elements progenitors, as an 
explanation of the large scatter, having been previously pro- 
posed by Ishimaru & Wanajo (1999) and Argast et al. (2000) 
(although with different mass ranges and mass dependences of 
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Table 2. The used prescriptions of the mass fractions of newly-produced elements for neutron capture elements in massive stars 
(r-process). 
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Fig. 2. [Eu/Fe] vs [Fe/H] . The abundances of simulated living 
stars at the present time are indicated by the blue dots. The red 
triangles are the observational data from various authors (see 
Table [B The dashed line is the prediction of the homogeneous 
model (Cescutti et al. 2006, model 1). 



i i i | i i i i | i i i i | i i i i | i i i r 



1 *A A A A 




[Fe/H] 

Fig. 3. As in Fig. 2, but for [Ba/Fe]. The authors of the obser- 
vational data are listed in Table \T\ The dashed line is the pre- 
diction of the homogeneous model (Cescutti et al. 2006, model 
1). 



the yields). Nevertheless, Ishimaru & Wanajo (1999) do not 
show results for ^-elements; Argast et al. (2000) show results 
for a-elements, but the resulting spread of their model is too 
large compared to the observed one. We have to emphasize 
that if we use, the original yields by Woosley & Weaver, rather 
than the yields by Francois et al.(2004), our model also shows 
a larger spread than the observed spread for the Mg (but not 
in the case of Ca). This is due to the strong dependence, in 
the original work of Woosley & Weaver, of the Mg yields on 
stellar mass, compared to the smaller dependence on the pro- 
genitor mass in the work of Frangois et al. (2004). However, 
we think that it is reasonable to use the yields of Frangois et 
al. (2004) for Mg, since recently many nuclear reaction data 
have changed, in particular in the case of this element. More 
recent nucleosynthesis yields, such as that found by Nomoto et 
al. (2006), show a smaller dependence on the progenitor mass 
as well. In the Figs. [TT] and [121 we show the relative frequency 
of stars at a given [El/Fe] ratio for different enrichment phases. 
To do that, we have used only the stars still alive in the halo, 
namely those with a mass < 0.8M o . The different enrichment 



phases: [Fe/H] < -2 and -2 < [Fe/H] < -1 are given in 
the panels from top to bottom for Si and Eu. In these Figs., the 
solid lines are the predictions of the model and the dashed lines 
are the observational data. 

At low metallicity, the model predicts a quite narrow dis- 
tribution for Si, whereas for Eu the distribution of the stars is 
broad, in agreement with the observations. At higher metallic- 
ities, both elements have predicted distributions slightly nar- 
rower than the observed ones. Nevertheless, we emphasize that 
the measured ratios are affected by observational errors, which 
in general are of the order 0.1 dex and could be responsible 
for the broader distributions of the observational data. In the 
low metallicity range, we noticed some discrepancies. For in- 
stance, the observed distribution of [Si/Fe] is slightly broader 
than the predicted distribution; therefore, a more careful anal- 
ysis of nucleosynthesis yields may be required (for instance, 
using metallicity dependent Si yields). The model is not able 
to reproduce ultra metal-poor stars, whereas stars with these 
metallicity are observed. In particular, the model predicts the 
formation of some metal free stars at the beginning, but then 
the metallicity jumps to [Fe/H] — 3.6. This may be due to an 
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Fig. 4. As in Fig. 2, but for [La/Fe]. The authors of the ob- 
servational data are listed in Table [T] The dashed line is the 
prediction of the homogeneous model (Cescutti et al. 2007a). 




[Fe/H] 

Fig. 6. As in Fig. 2, but for [Y/Fe]. The authors of the observa- 
tional data are listed in Table [T] The dashed line is the predic- 
tion of the homogeneous model (Cescutti 2007b). 



overly rapid rise of [Fe/H], related to the adopted infall rate. 
Indeed, numerical simulations (Recchi et al. 2001) show that 
the mixing timescale of the ejecta, although fast, is of the order 



of lOMyr, therefore larger than the assumed timestep. If in- 
cluded, this delayed mixing might increase the amount of very 
metal-poor stars. Moreover, the model predicts that ~ 5% of 
the simulated stars are metal free. This is the fraction of stars 
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Fig. 5. As in Fig. 2, but for [Sr/Fe]. The authors of the observa- 
tional data are listed in Table [T] The dashed line is the predic- 
tion of the homogeneous model (Cescutti 2007b). 
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Fig. 7. As in Fig. 2, but for [Zr/Fe]. The authors of the observa- 
tional data are listed in Table [T] The dashed line is the predic- 
tion of the homogeneous model (Cescutti 2007b). 
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Fig. 8. As in Fig. 2, but for [Si/Fe]. The authors of the observa- 
tional data are listed in Table [T] The dashed line is the predic- 
tion of the homogeneous model (Francois et al. 2004). 
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Fig. 10. As in Fig. 2, but for [Mg/Fe]. The authors of the ob- 
servational data are listed in Table [T] The dashed line is the 
prediction of the homogeneous model (Francois et al. 2004). 



formed in the first 5 Myr, in which no SN has yet exploded and 
enriched the ISM. We emphasize that the star formation can be 
slowed down, for example, by assuming a different infall rate 
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Fig. 9. As in Fig. 2, but for [Ca/Fe]. The authors of the ob- 
servational data are listed in Table Q] The dashed line is the 
prediction of the homogeneous model (Francois et al. 2004). 



or a strong outflow of gas from the halo, as recently adopted by 
Goswami & Prantzos (2000), which would modulate a different 
star formation. This can solve the absence of ultra metal-poor 
stars and the too large fraction of metal free stars. The problem 
of the metal free stars can be also solved by adopting a dif- 
ferent IMF. The results of many works (see Larson 1998; Abel, 
Bryan, & Norman 2000; Hernandez & Ferrara 2001 ; Nakamura 
& Umemura 2001; and Mackey, Bromm, & Hernquist 2003) 
predict that stars formed in a metal free gas must be massive 
stars. The long-living stars, which are low mass stars, start to 
form when the ISM has been already enriched by these massive 
and metal free stars, the so-called Population III. The effect of 
the Pop III on the global chemical evolution should be negligi- 
ble (see Matteucci & Pipino 2005, Matteucci & Calura 2005), 
as it involves a small total amount of recycled mass. In fact, we 
tested a model with a top-heavy IMF, for the very first period, 
up to a metallicity Z = 10" 4 Z o ; we decided to use this value 
because in the literature 10" 4 Z o is considered the metallicity of 
the transition to have a normal IMF (see Schneider et al. 2002). 
This stage lasts only about 5-10Myr, depending on the chosen 
cubic region. This top-heavy IMF has the same shape of the 
one we normally use (Scalo et al. 1986), but we changed the 
lower limit from 0.1 M Q to 1M Q . We found, in this way, that the 
chemical evolution does not change but we eliminated the pres- 
ence of metal-free stars, which are not observed. In fact all the 
stars, born during this stage will die before the present time and 
the stars born after will not be metal-free. We emphasize that 
we have decided to use the same parameters and conditions of 
the homogeneous model to test the validity of the inhomoge- 
neous one. It is clear from the Figs, that this new model for a 
high number of star formation events well approximates the ho- 



8 



G. Cescutti: Inhomogeneous model for the Galactic halo 



KJ.2 



1 1 1 1 1 1 1 

—J 

1 1 1 1 r— | ri 


i i i i i 
[Fe/H] < -2 J 


r 

: ! 


— 1 1"!-^^— i 1 1 1 

■2 < [Fe/H] < -1 J 
"In : 

■n h - 



[Si/Fe] 

Fig. 11. The relative frequency of stars at a given [Si/Fe] ratio 
for different enrichment phases. In solid line the predictions of 
the model for the living stars at the present time, in dashed line 
the observational data. The authors of the observational data 
are listed in Table [T] 



mogeneous model, which is shown by the dashed line in Figs. 
2-10. 



5.2. The ratios [Ba/Eu] and [Ba/Y] 

Our model predicts a spread in the abundance ratios of two 
elements if the sites of production of these two elements are 
different, or if the ratio of the yields of these same elements 
presents large variations as a function of the stellar mass. 

This is not the case for the yields of neutron capture ele- 
ments. In Fig. [13] we show the ratio [Ba/Eu] versus [Fe/H] . As 
expected, the results of the model shows a small spread, as for 
the plot of [or/Fe], being the ratio between the newly-produced 
Ba and the newly-produced Eu almost constant as a function 
of the stellar mass. On the other hand, the observational data 
show a spread larger than predicted in the metallicity range 
-2 < [Fe/H] < -1. We observe that most of the data seem 
to have a larger overabundance of Ba than that predicted by our 
model. Taking into account this fact, the observational spread 
could be explained either by an earlier production of B a by s- 
process in intermediate stars from 3 to 8M© (that we do not 
take in account in our nucleosynthesis); by self enrichment of 
the observed stars due to dredge-up; or, to a binary system with 
mass transfer from an AGB star, which is not shining anymore, 
to the presently observed companion star (see Aoki et al. 2006). 

In Fig. [14] we show the ratio [Ba/Y] versus [Fe/H]. The 
results of the model relative to this ratio are not satisfactory. 
The observational data show a very large scatter at [Fe/H] ~ 
-3 that is not predicted by our model. Contrary to the [Ba/Eu] 
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Fig. 12. As in Fig. EH but for [Eu/Fe]. The authors of the ob- 
servational data are listed in Table [T] 



ratio, for which the available data show a moderate spread at 
[Fe/H] ~ -3, [Ba/Y] ratio shows a large spread. One possible 
explanation for this spread could be that the r-process yields 
that we use in our model are indicative of the mean contri- 
bution by r-process to the abundances of these elements. We 
emphasize that these two elements are in different peaks in the 
solar abundances, both for what concern s-process (not so im- 
portant at this stage) and for what concern r-process. It could be 
that, as introduced by Otsuki et al. (2003), the r-process is not 
unique but consists of different contributions and what we use 
are only the mean values. Massive stars may produce r-process 
with different patterns, probably as a function of the multiple 
factors that influence r-process. In our inhomogeneous model, 
we use only one pattern for all the neutron capture elements. 
This could be the reason why we are not able to reproduce the 
spread for the ratios of neutron capture elements. Moreover, 
this problem should be more visible when we compare neutron 
capture elements belonging to different peaks as in the case of 
[Ba/Y]. 

6. Conclusions 

In this paper, we tried to solve the problem of the large spread 
observed at low metallicity of the [el/Fe] ratios for neutron cap- 
ture elements and, at the same time, the smaller spread ob- 
served in [or/Fe] ratios. We developed a new model for the 
chemical evolution of the halo of the Milky Way. In this model, 
we coupled the birth of stars with random masses, with dif- 
ferent mass ranges for the production of or-elements and neu- 
tron capture elements. In particular, the site of production of 
or-elements is the whole range of the massive stars, from 10 to 
80M© whereas the mass range of production for neutron cap- 
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Fig. 13. [Ba/Eu] vs [Fe/H]. The abundances of simulated liv- 
ing stars at the present time in blue dots, the red triangles 
are the observational data from Burris et al.(2000); Francois 
et al. (2007); Fulbright (2000, 2002); Gilroy et al. (1998); 
Gratton & Sneden (1994); Honda et al. (2004); Ishimaru et al. 
(2004); Johnson (2002); Mashonkina & Gehren (2000, 2001); 
McWilliam et al. (1995); McWilliam & Rich (1994); Prochaska 
et al. (2000); Ryan et al. (1991, 1996); Stephens (1999); and 
Stephens & Boesgaard (2002). The dashed line is the predic- 
tion of the homogeneous model (Cescutti et al. 2006, model 
1). 



ture elements lies between 12 and 30M o . We showed that these 
two assumptions can explain the larger spread in the abun- 
dances of metal-poor stars for neutron capture elements and the 
smaller spread for a-elements. We note that the idea of a small 
range of the r-elements progenitors, as explanation of the large 
scatter, is not new, having been already proposed by Ishimaru 
& Wanajo (1999) and Argast et al. (2000), although the mass 
ranges and the mass dependences of their yields and this work 
are different. We emphasized that the parameters we used (de- 
scribed in Sect. are the same as the homogeneous model (see 
Cescutti et al. 2006). We adopted this point of view because 
these parameters have been already constrained to give good 
results compared to the observational data at higher metallic- 
ities and compared to the mean trends. In fact, toward high 
metallicities ([Fe/H]>-2.0 dex) the model naturally gives re- 
sults compatible to the homogeneous model. However, this set 
of parameters is a starting point and the model still needs a 
better investigation of the parameter space. In particular, con- 
cerning the early galactic stages, the model is not able to repro- 
duce ultra metal-poor stars, even though stars with this metal- 
licity are observed. This may be due to an overly rapid rise of 
[Fe/H], related to the used infall law. To avoid these problems, 
a different infall law can be used or an outflow from the halo 
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Fig. 14. As in Fig.Obut for [Ba/Y] vs [Fe/H]. The data for this 
elements are from Burris et al.(2000); Francois et al. (2007); 
Fulbright (2000, 2002); Gilroy et al. (1998); Gratton & Sneden 
(1994); Honda et al. (2004); Ishimaru et al. (2004); Johnson 
(2002); McWilliam et al. (1995); McWilliam & Rich (1994); 
Nissen & Schuster (1997); Prochaska et al. (2000); and Ryan 
et al. (1991, 1996). The dashed line is the prediction of the 
homogeneous model (Cescutti 2007b). 



can be included to provide a smoother increase of the metal- 
licity. Moreover, this new model generates an over abundance 
of metal free stars if we do not impose at least a slightly dif- 
ferent IMF for the first episodes of stars formation to solve the 
problem of the metal free stars. 
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